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ABSTRACT 

Magnetar flares generate Alfven waves bouncing in the closed magnetosphere with energy up to 
~ 10^6 erg. We show that on a 10-ms timescale the waves are transmitted into the star and form a 
compressed packet of high energy density. This packet strongly shears the stellar crust and initiates 
a plastic flow, heating the crust and melting it hundreds of meters below the surface. A fraction of 
the deposited plastic heat is eventually conducted to the stellar surface, contributing to the surface 
afterglow months to years after the flare. A large fraction of heat is lost to neutrino emission or 
conducted into the core of the neutron star. 

Subject headings: dense matter — magnetic fields — stars: magnetars — stars: neutron — waves 


1. INTRODUCTION 

Magnetars are luminous slowly rotating neutron stars 
powered by the decay of ulstrastrong magnetic field s 
B = 10^"^ — lO^^G (see e.g. IWoods &: Thompsonl (1200611 : 
iMereghettil (1200811 for reviews). They have hot sur¬ 
faces, produce nonthermal magnetospheric radiation and 
strong bursts of hard X-rays. Occasionally, magnetars 
produce giant flares with energies of 10“^^ — 10^® ergs. To 
date three giant flares have been observed from three 
magnetars. The main peak of the giant flare lasts ~ 0.3 s 
and can reach huge luminosities ^ ~ 10^^ erg s“^. Less 
powerful flares with < 10^^ erg s“^ (often called 
“bursts”) occur much more frequently. 

The flares are associated with a sudden change in 
the magnetospheric configuration, which could be trig- 
gered by an instability inside or outside the neutron star 
(|ThomDSon fc DuncanlllOOOl . This cataclysmic event in¬ 
volves strong deviations from the magnetostatic equi¬ 
librium, launching waves of large amplitudes. Part of 
the released magnetic energy is promptly dissipated and 
converted to radiation, and part is stored in the excited 
waves. 

In particular, Alfven waves are generated with a to¬ 
tal energy up to ^ 10“^® erg. They are trapped on 
the closed magnetic field lines, as the group velocity of 
Alfven waves is parallel to the magnetic field. The fate 
of their energy is poorly known. It was proposed that 
the Al fven waves can be damped through nonlinear pro¬ 
cesses (|ThomDSon h BlaeillQQSH , which become efficient 
at very large amplitudes of the waves (see Section 5.2). 

In this paper, we propose another mechanism of the 
Alfven wave dissipation, which results from the wave in¬ 
teraction with the star. The waves are ducted along the 
magnetic field lines with nearly speed of light and reach 
the stellar surface on a millisecond timescale. In Sec¬ 
tion 2 we examine the wave interaction with the star 
and find that a significant fraction of the wave energy is 
transmitted into the stellar crust. The reflected waves 
keep bouncing in the magnetosphere, however in a few 
tens of milliseconds most of their energy is drained and 
deposited into the crust, in the form of a compressed 
shear wave packet. In Section 3, we show that this packet 


causes strong plastic heating of the crust. In Section 4 
we investigate the fate of heat deposited by the plastic 
damping of Alfven waves. In particular, we evaluate the 
heat flux conducted back to the surface and the result¬ 
ing surface luminosity, which should emerge long after 
the flare. Our results are discussed in Section 5. 

2. WAVE TRANSMISSION INTO THE GRUST 

The crust is nearly incompressible and supports shear 
waves which can be excited by the Alfven waves im¬ 
pinging from the magnetosphere. Excitation of two-fluid 
crustal modes can be neglected, and the recent claim 
that magnet ospheric Alfven w aves transform into crustal 
Hall waves (iLvutikovI I20I5I1 is incorrect. Hall waves 
propagating parallel to the magnetic field B with fre¬ 
quency u! in the crust of density p have refraction index 
Af = ck/uj = ujpej ~ ^^(we use 

the standard notation X^, = V/IO'" for a quantity X in 
cgs units). This implies a huge impedance mismatch with 
the magnetospheric Alfven waves, which have Af Ri 1, 
and therefore their transformation to Hall waves is sup¬ 
pressed Q No significant separation between electron and 
ion velocities can occur on ms timescales, and the re¬ 
sponse of the crust to the external disturbance is essen¬ 
tially single-fluid. 

2.1. Transmission coefficient 

Gonsider a magnetospheric Alfven wave of frequency uj 
impinging on the crust of the neutron star. For simplic¬ 
ity let us assume that the initial (unperturbed) magnetic 
field Bz is uniform and vertical, so the wave is prop¬ 
agating vertically along the z-axis, and the horizontal 
displacement ^( 2 ) is along the y-axis. The plasma-filled 
magnetosphere and the crust are excellent conductors; 
therefore the magnetic field is frozen in the medium and 
the horizontal field By is related to the displacement by 

^ Only electrons move in a Hall wave (analogous to whistler in 
plasma physics) while ions are static. The velocity of the electron 
fluid vh = j/eue is related to electric current j = (c/47r)V X B, 
which gives a tiny uh because of the high electron density rie in the 
crust. Therefore, the “two-fluid” (electron-ion) description is useful 
only for slow phenomena in the crust. 
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BylBz = d^jdz. The wave speed in the magnetosphere 
is close to the speed of light c, and the wavelength is 
Ao = 27rc/w. 

The wave propagation is described by the equation, 



Bl 

47rc^ 




Bl ^ 

47r 9z2 ' 


( 1 ) 


Here p{z) is the mass density and p + BHAttc? can be 
thought of as the effective inertial mass density of the 
magnetized medium^ The first term on the right-hand- 
side describes the restoring force of magnetic tension, and 
the last term describes the force due to the shear stress 
in the medium. 

In particular, if the medium is elastic with a shear 
modulus p then a = —ps, where s = d^fdz is the strain 
of the elastic deformation. In this case. Equation m 
becomes a simple wave equation with the wave speed 
given by 


, ^ B^/Att + p{z) 

^ ’ BI/Attc^ +p{zy 


( 2 ) 


In the magnetosphere, we will neglect the mass den¬ 
sity p and the shear modulus p, which gives v = c. 
In the crust, we will use the profiles p[z) and p{z) 
shown in Figure [TJ The density profile is obtained from 
the relati vistic hydrostatic equation using SLy equation 
of state (jHaensel fc Potekhinl 120041 1 for a neutron star 
with mass M = I. 4 M 0 I 3 The radius of the star is 
R = 11.7 km, and its surface gravitational acceleration 
is g = {GM/R^){1 — rgjR)~^^'^ = 1.7 x 10^^ cm s“^ 
where rg = 2GM/c^. For the shear m odu lus u we use 
the ht ting formula given bv iPir r] (l200l and lSotani et all 
(j2007tl for low and high densities. 

As the wave propagates into the deeper crust, its speed 
is reduced and its wavelength is compressed. 


A(z) = Ao 
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( 3 ) 


The reflection of the wave occurs in the region where 
the characteristic scale-height for the change of v{z), 


h{z) 


V 

\dv/dz\ ’ 


( 4 ) 


is smaller than the wavelength A(z). Figure [5] shows A(z), 
h{z), and the depth zi where they are equal. The typical 
value of Zi is around 200 meters below the surface; its 
exact value depends on B^. 

The transmitted wave below Zi has X h and can be 
described in the WKB approxim ation. Then the w ave 
displacement takes the form (e.g. iFitzoatridi (I2013IF . 




const 

Zl/2(^) 
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( 5 ) 


^ This expression is approximate as it neglects the contribution 
from the horizontal field component By. In the models presented 
below, By > Bz when the wave propagates into the crust; however, 


this only occurs in the dense region where B^ jAn pc^ and the 
magnetic field inertia anyway may be neglected. In the region 
where B^ f An ^ pc^ the wave amplitude s = ByjBz < 1 and it is 
acceptable to approximate B^ B^. 

^ The ultrastrong magnetic field significantly changes pressure 
where the electron Fermi energy is below the Landau energy. This 
impacts the density profile p ( z ) at shallow depths. However, at 
depths of interest in this paper (where p ^ 10^ g cm“^) this effect 
is small and neglected. 



Fig. 1. — Density p { z ) and shear modulus pi { z ) of the neutron 
star crust. The star is assumed to have mass M = I.AMq and SLy 
equation of state. 



Fig. 2. — Solid curve shows the wavelength of the shear wave 
propagating in the magnetized crust, A = Xqv/c, where Aq = 10 km 
and V is the speed of the wave. Dashed curve shows the character¬ 
istic scale-height of the wave deceleration, h = v/\dzv\‘ A vertical 
magnetic field Bz = 3 x 10^^ G is assumed in this example. 


where Z{z) is the impedance 
Z{z) = 


\ B'i , J 

1/2 

[52 , 1 



[4^02 
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( 6 ) 


A simple estimate for the trans mission coefficieiit is 
obtained using the impedance at Zi (IBlaes et al.lll989ll . 


^ 4Z{zyZ{0) _ 4^ 

^ [Z(zi) + Z(0)f ^ C ■ 

For instance, for = 3 x 10^^ G, Equation ([7]) gives 
T ^ 5%. A more accurate transmission coefficient is ob¬ 
tained by solving numerically the wave equation, which 
gives a higher value of T = 12% (the smoothness of the 
crustal density variation as the wave approaches zi en¬ 
hances the transmission). The numerically calculated 
T{Bz) is shown in Figure |31 It is comparable to O.I for 
typical magnetar fields. 

The reflection coefficient 77. = I — T 0.9 is large, and 
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Fig. 3.— Transmission coefficient 7” as a function of vertical 
magnetic field Bz , obtained from the numerical solutions of Equa¬ 
tion ifTIl . 


the reflected Alfven waves will bounce many times in the 
magnetosphere. Their amplitudes decrease by T ~ 10% 
every time they bounce from the surface. The repeated 
transmission events form a train of compressed waves 
in the crust. This train propagates into the crust with 
velocity v ~ 10“^c. 

One can show from Equation ([5]) that the strain 
s = d^/dz in the transmitted wave evolves as |s| oc 
oc It increases as the wave propagates 

into the deeper and denser crust. This has a simple phys¬ 
ical reason: the wave decelerates, and hence its energy 
density grows as v~^ (so that the wave continues to 
carry its energy flux = U^v = const). The wave 
energy oscillates between the kinetic energy and the hor¬ 
izontal field plus elastic energy of the crust. Therefore, 
may be written in two ways: Uw ^ {p+B'^/A ttc^) 
(kinetic) or Uw ~ s"^ /An + p) (magnetic-|-elastic). In 

the region where pc^ > B'^/An and p < B^/An this re¬ 
quires oc (you)”^ and s^ oc v~^. Then the relation 
s~^/Aoc^/u gives 

^ oc V oc s oc (8) 

In the lower crust where p ^ B^/An and v k. i^pjp)^!"^ Ri 
const m 10® cm s“^ one finds s oc p~^/‘^. In this region 
the wave strain significantly decreases. This evolution of 
s with depth (increase and then decrease) may be ob¬ 
served in the numerical simulation presented below. 

2.2. Numerical model 

To illustrate the transmission process we set up a sim¬ 
ple one-dimensional simulation of waves bouncing in the 
magnetosphere between the footprints of a closed mag¬ 
netic flux tube. The Alfven waves are ducted along the 
magnetic field lines and the problem can be made one¬ 
dimensional by pretending that the flux tube is straight, 
and by placing its two opposite footprints on the z-axis, 
separated by distance L. Here L represents the length 
of magnetospheric field lines. The stellar crust with 
the density profile p{z) is placed symmetrically at the 
two ends of the computational box. The crust thickness 
^ 1 km is much smaller than L. 


Fig. 4.— A snapshot of the wave after four refiec- 
tion/transmission events. Solid red curves show the magnetic field 
lines deformed by the horizontal displacements in the wave; back¬ 
ground grey color shows the density of the crust. The wave front 
has reached the depth of 800 m by this time. One can see the four 
oscillations in the transmitted and compressed packet. The shape 
of each oscillation reflects the initial shape of the magnetospheric 
wave assumed in the simulation. This simulation included no plas¬ 
tic dissipation and assumed the magnetic field Bz = 3 x 10^“^ G. 


For any initial shear distortion of the field lines, one 
can calculate the subsequent dynamics of the generated 
Alfven waves by solving numerically Equation ([T]). In 
our numerical models, we take the initial distortion of 
the form, 


^o(z) = A exp 


{z - zo)^ 
2P 


( 9 ) 


which is localized in the middle of the box zq, far away 
from the crust, with I < L. The distortion immedi¬ 
ately splits into two waves propagating toward the op¬ 
posite ends of the box. The strain profile of each wave 
s{z) = d^/dz is determined by the initial distortion. An 
important parameter of the wave is its initial maximum 
strain. 

So = i max|a^^o(z)| = (10) 

The total energy initially stored in the two waves is 




( 11 ) 


where S is the cross section area of the flux tube. 

We follow the evolution of the waves and their inter¬ 
action with the crust until almost all wave energy Eq 
has been drained from the magnetosphere; this typically 
takes tens of light-crossing times L/c. The wave equa¬ 
tion o is solved on a grid with 1000 points in the mag¬ 
netosphere (uniformly spaced) and a much finer grid in 
the crust (one point per meter). Convergence tests have 
been done to ensure that the grid is sufficiently large to 
resolve the wave dynamics. 

The snapshot of the simulation in Figured shows the 
distortion of the crust at time t = A.6L/c, when the mag¬ 
netospheric waves have bounced four times. The param¬ 
eters of this sample model are L = 40 km, I = 5/-\/2 km, 
A = 5 km, and 5^ = 3 x 10^^ G. The corresponding 
So = (2e)“^/^ Ri 0.43 and Eq/S A.5 x 10®^ erg cm“^. 
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In the snapshot shown in the figure, about 1/3 of the 
wave energy Eq has already been transmitted into the 
crust. The transmitted wave has been decelerated to 
V Ri 10® cm s“^ and compressed by the factor of c/v k. 
3 X 10^. The compression creates a high energy den¬ 
sity of the horizontal magnetic field aX z = 200 — 500 m, 
{sBzY/S-K ^ 10(soi3z)^/87r, and strain s ^ 3so- 

3. PLASTIC HEATING 

The description of the wave dynamics in Section [5] 
is incomplete because it assumes the elastic response 
a = —/rs everywhere in the crust. The more realistic 
model must take into account two facts: (1) When the 
solid crust is deformed by the shear wave beyond a crit¬ 
ical stress O’er its response becomes plastic rather than 
elastic. (2) The crustal temperature may be high enough 
to reduce Ocr or even melt the crust, leading to Ocr ~ 0. 
Therefore, the model should keep track of the crustal 
temperature. 

3.1. Pre-flare temperature profile 

The typical persistent sur face temperature of magne - 
tars is Tg ^ (3 — 4) x 10® K (|Woods fc Thompsonll2006ll . 
It corresponds to the radiation energy flux F = aT^ ^ 
10^^ erg s“^, where a = 5.67 x 10“® erg s“^ cm“^ is 
the Stefan-Boltzmann constant. A usual way to estimate 
the subsurface temperature profile of neutron stars T(z) 
assumes that the surface flux F is supplied by quasi¬ 
steady diffusion of heat from the crust. Then T(z) is 
given by the equation, 

«(T,z)—= F = o^,^ (12) 

where k is the effective conductivity, which is dominated 
by degenerate electrons at densities p > 10® g cm“® 
and by radiation in the low density layers near the sur¬ 
face. Note that k depends on the local magnetic field. 
Combining Equation (TT^ with the hydrostatic equation 
dP/dz = pg gives 

dlogT 3 PKT^ 
dlogP 16 5 

Here P is the pressure, g = [GM/R^){1 — TglR)~^^‘^ 
is the surface gravitational acceleration, and K = 
/Hup is the effective opacity; all quantities are mea¬ 
sured in the local rest frame of the crust. The surface 
luminosity and temperature measured by a distant ob- 
server are .if^ = (l-rg/i?).if and = {l-rg/Rf/'^T, 
(lThornelll97^ . 

Equation m assumes that the temperature profile 
had enough time to relax to the steady state at depths of 
interest, which typically takes ^ 1 yr. In a true steady 
state, the relatively high surface temperature of magne- 
tars requires a source of heat at d epths of a few hun¬ 
dred meters (iKaminker et al.l 1200611 . Alternatively, one 
may view this temperature profile as qausi-steady, slowly 
cooling after a previous heating episode. 

If one accepts this thermal model for the pre-flare state 
of the crust, one can find T{z) from Equation (fOl) and 
determine the melting depth Zm sit, above which t he crust 
is melted. We use the code of iPotekhinI (|1999ll to cal¬ 
culate the thermal conductivity and the melting point 


Tmeitip) of the crustal material. An approximate result 
is sufficient for the purposes of this paper and we do not 
discuss here the poorly known chemical composition of 
the magnetar crust. For simplicity, we assume an iron 
crust with a small impurity parameter. 

We also assume that the pre-flare magnetic field is 
not far from vertical. This is a reasonable assumption 
for the melted layer {z ^ 100 m, see below) where the 
crustal magnetic field should match the force-free mag¬ 
netosphere. A strong toroidal field can only be stored in 
the core of the neutron star or the deep crust. 

Equation (fT^ can be solved numerically as described 
in detail in previous works, which calculated the rela¬ 
tion between the surface effective temperature Tg and 
the internal temperature Th measured at neutron-drip 
dept h Zh zz 400 m (ph = 4 x 10^ ^ g cm“®). In partic¬ 
ular, IPotekhin fc YakovlevI (1200111 provide a fitting for¬ 
mula for Tb(Ts) for various magnetic fields and assuming 
an iron crust. We use this relation to impose the condi¬ 
tion T = Ti, aX z = Zb, and then reconstruct the profile 
T{z) in the region of interest p > 10® g cm“® (above or 
below Zb) by integrating Equation (fl^ from Zb- Thus 
we avoid integration in the shallow surface layers where 
thermal conductivity is dominated by radiation, and so 
we only use the electron conductivity in our calculations. 

For a given Tg, this calculation gives the subsurface 
temperature profile T{z) and the melting depth Zmeit- 
The melting temperature is approximately given by 

T^eit ~ 2.4 X lOVlf K. (14) 

The exact value of melting depth Zmeit depends on the 
magnetic field and its orientation relative to the stellar 
surface. A strong field increases the thermal conductivity 
along B and decreases it perpendicular to B. Therefore, 
a vertical field tends to reduce the internal temperature, 
thus decreasing Zmeit- A horizontal field would hamper 
the heat flow in the vertical direction and increase Zmeit- 
A typical Zmeit in magnetars with non-horizontal surface 
fields is ^ 100 m. 


3.2. Plastic flow 

As the wave packet propagates into the crust below 
2^meit, it starts to interact with the solid phase (lattice). 
The response of the lattice is elastic as long as its strain 
is below a critical value Scr- The maximum Scr ~ 0.1 is 
comparable to the yielding th reshold for an ideal crystal 
(|Chugunov fc Horowitd 1201011 . The actual strain in the 
wave s ~ so(7'c/v)^^‘^ is much higher than Scr, and so 
the wave initiates a strong plastic flow with the high 
frequency w. In contrast to fluid Alfven wave or elastic 
shear wave, the plastic flow is dissipative, i.e. it converts 
the wave energy to heat, reducing its amplitude. Below 
we include this process in our wave propagation model. 

The plastic heating rate per unit volume is 


dt 




(15) 


where Spi = s — Sei is the plastic part of the strain, Sei 
is the elastic part, and a is the shear stress sustained by 
the plastic flow. In the plastic regime \a\ > a^r where 
(Jer = f i Srr- T he simple model of “viscoplastic solid” (e.g. 
Ilrgensl (1200811 1 gives the stress of the plastic flow in the 
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form, 


|cr| = CTcr + ?7|Spl|, (16) 


where 77 is a viscosity coefficient. 

The crystal becomes “soft” (i.e. dcr drops) if it is 
heated to a temperature comparable to the melting point 
Tmeit- The softening effect is responsible for the ther¬ 
moplastic instability that can release internal ma gnetic 
stresses in magnetars (|Beloborodov fc Levinl[20l^ . This 
instability however develops on a timescale much longer 
than 10 ms and does not affect the dynamics considered 
in this paper. Here the plastic flow is driven by the strong 
external magnetic stress from the flare (rather than de¬ 
velops spontaneously inside the crust) and immediately 
reaches huge strains |s| ^ Scr and high temperatures. 
Because |s| ^ Scr, the detailed behavior of Scr(T) and 
o'cr(T) is not important; our calculation should merely 
take into account the fact that plastic heating switches 
off when T approaches Tmeit- 

This effect is included as follows: the stress a of the 
plastic flow is multiplied by the factor 1 — C/th/dmeit, 
where Uth{p,T) is the thermal energy density and 
Clmeit = C4h(p, Tmeit)- This prescription enforces cr = 0 
when T = Tmeit- 

The stress in the elastic regime a = —fis must match 
the plastic stress at s = Scr- This condition is automati¬ 
cally satisfied for the cold crystal. For a hot crystal the 
reduction of (Tcr(T) may be interpreted as the reduction 
of shear modulus p or the reduction of Scr (or both). 
These details are not important for our model, because 
the plastic flow has |s| ^ Scr- The numerical models pre¬ 
sented below assume Scr(T) = 0.1 = const and use the 
following prescription for a, 


( t^th \ f elastic 

CT = ( 1- I X < (0.1/r -I- ? 7 Spi) sign(—s), plastic 

V Ume\t) y 0, liquid 

(17) 

where p is the shear modulus at T <g; Tmeit shown in Fig¬ 
ure 1 , and one may think of (1 — Ut\i/Unie\t)p as the shear 
modulus reduced by heating. We verified that practically 
the same results are obtained if we choose a temperature- 
dependent Scr = 0 . 1(1 — C/th/C7meit) with shear modulus 
unchanged by heating. 

Finally, we must choose 77 , which is unknown for the 
crustal material. The transition between the plastic and 
elastic regimes is smooth if 77 vanishes when |s| = Scr. 
Therefore, we assume 77 of the form 77 = a/r | |s| — Scr|, 
where a is a constant. We tried various values of a and 
found that plastic heating weakly depends on it as long 
as a is sufficiently large, a > 3 x 10“® s. Our sample 
numerical models use a = 3 x 10“^ s. 

The dynamic system described by Equation ([T]) with a 
given by Equation dig satisfies the energy conservation 
law. 


d? “ ^ / dt “ dt + Eb+E,i), (18) 



Fig. 5. — Shear wave propagation in the magnetar crust viewed 
on the spacetime diagram: Upper panel: horizontal displacement 
of the wave, Middle panel: strain s = d^/dz. Lower panel shows 
where the crust is deformed elastically (white), flowing plastically 
(black), and melted (gray). 

where 


i^kin = ^ J 

'(-.so 


(19) 
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where |sei| < Scr in the elastic zone and |sei| = Scr in the 
plastic zone. 

The plastic flow occurs where \a\ exceeds and con¬ 
tinues as long as d\s\/dt > 0. Whenever the local ab¬ 
solute value of the strain stops growing, the plastic flow 
switches to the elastic regime; at this point Sei and a are 
reset to zero. 

3.3. Wave damping and post-flare crustal temperature 

We re-run the model described in Section 2.2 with the 
new expression for a that takes into account the plastic 
damping in the crust (Equation (flTl) '). The initial state is 
assumed to have the surface temperature Tg = 3 x 10® K. 
All other parameters are the same as in Section 2.2, in 
particular = 3 x G and sq = (2e)“^/^ ~ 0.43. 
The spacetime diagram of the wave evolution in the crust 
is presented in Figured It shows the wave displacement 
and strain, and indicates the elastic, plastic, and melted 
regions. 

Figureinishows the history of the wave energy transmis¬ 
sion from the magnetosphere to the crust and the plastic 
damping effect. One can see that most of the transmit¬ 
ted wave energy is promptly converted to heat. We have 
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Fig. 6.— Solid curve shows the evolution of the energy fraction 
left in the magnetosphere and dashed curve shows the energy frac¬ 
tion converted to heat. Each step in the solid curve corresponds 
to the simultaneous reflection of the two symmetric waves bounc¬ 
ing in the magnetosphere in the opposite directions. The reflec¬ 
tion coefficient 'K. 0.88, and the magnetosphere loses energy as 

where N is the number of reflection events. Each step takes 
L/c = (4/30) ms, so Bmag/Eo ^ 0.88*/°133ms 


verified that our numerical simulation satisfies the con¬ 
servation law (Equation (TTSl) ') with accuracy better than 
1%. Each time the wave hits the surface, the transmis¬ 
sion coefficient is approximately 12%, and almost all the 
wave energy Eq is damped after ^10 ms, so most of Eq 
becomes stored as crustal heat. This heating results in 
deep melting of the crust, down to 500 m. 

Since plastic dissipation switches off at the meting 
point, the crust naturally acquires the “ceiling” temper¬ 
ature T Ri Tm in an extended region below the surface. 
The resulting temperature profile immediately after the 
flare is shown in Eigure[7l To investigate how the results 
depend on and sq; we have calculated the models 
with G = 0.3,1,3,10 and sq = 0.13,0.25,0.43. 

Stronger waves in stronger magnetic fields melt deeper 
layers of the crust, up to 600 m in the calculated models. 


4. COOLING 

After the flare, the hot crust will cool on a much longer 
timescale. Two main processes cool the crust: neutrino 
emission and heat conduction. The temperature evolu¬ 
tion with time is described by the following equation. 



^ dz J 


- qy, 


( 22 ) 


where k, is the thermal conductivity and Cv is the heat 
capacity of the crust; both are functions of local p{z), 
T(z,t), and B. The sample numerical models presented 
below assume a uniform vertical magnetic field B = = 

const. We u se k(p,T B) and Cv{p,T, B) calculated by 
the code of iPotekhii] (|1999l) . The term q^{p,T,B) is 
the rate of local cooling b y neutrino emission. This rat e 
is described in detail by IPotekhin fc Yakovlev! (|2001l) . 
They provide useful analytical approximations for four 
relevant channels of neutrino emission: plasmon decay, 
bremsstrahlung, synchrotron, and electron-positron an¬ 
nihilation. We use their formulas in our calculations. 



Fig. 7.— Temperature profile after the magnetospheric Alfven 
waves have been absorbed by the crust; upper panel for fixed sp = 
0.43 and different Bz, and lower panel for fixed Bz = 3 x G 
and different so- The left boundary in the figure is chosen at depth 
2 : 60 m where p = 10^ g cm“^; at 2 ^ 100 m the crust is melted 

before the flare, and hence no plastic heating can take place. 


Similar to previous simulati ons of time-dep e ndent 
heat diffusion in magnetars ([ Kamink er et al.l 120061 
iBrown fc Cummin5l20o4 iPons et al.ll2009[) .~ we separate 
the crust into two regions: a blanketing envelope and an 
interior region. Here we choose the envelope boundary 
at Zfc Ri 60 m where p = pb = 10® g/cm^. The typical 
timescale of heat diffusion from this depth is tt 10®s. 
It is sufficiently short to give a quasi-steady state in the 
envelope, and so the steady-state solution may be used 
to determine the relation between Tb = T{zb) and the ef¬ 
fective surface temperature T^. Note that Tg defines the 
energy flux E = oTf through the envelope z < Zb, and 
thus in essence the Tb-Tg relation is a relation between 
Tb and the heat flux F = ndT/dz at Zb- It serves as 
a boundary condition for our time-dependent heat diffu¬ 
sion problem &t z > Zb- Since this boundary condition 
relies on the steady-state solution z < Zb it can only be 
accurate when T{t,z > Zb) evolves on timescales longer 
than tb- 

We calculated the Tb-Tg relation at pb = 10® g cm“® 
using the steady-state solutions obtained in Section 3.1. 
This gave a tabulated boundary condition F{Tb) at the 
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Fig. 8.— Evolution of the crustal temperature profiles in our 
fiducial flare model with B = Bz = 3 x 10^'^ G and sq = 0.43. 

upper boundary of our computational box Zb ~ 60 m. 
The lower boundary is chosen at z ~ 1 km, near the 
bottom of the crust where p ~ 10^^ g cm“^. The exact 
position of the lower boundary is not important as long 
as it is deep enough. The deep crust has a high thermal 
conductivity and the heat is absorbed by the (approxi¬ 
mately isothermal) core of a huge heat capacity. We use 
the absorbing boundary condition at constant tempera¬ 
ture T ^ 3 X 10® K, neglecting the increase of the core 
temperature due to the absorbed heat. 

The initial condition T^z^Q) for Equation (I22|) is pro¬ 
vided by the plastic heating model (Figure [7|). The ini¬ 
tial temperature increases along the melting curve Tm{z), 
reaches maximum, and drops at larger depths. We evolve 
this initial temperature profile on a uniform grid of 1000 
points and a time step of 10 s for 10^ steps. To speed up 
the simulation, the values of k, Cy, and on the grid 
are updated every 500 time steps. Convergence tests ver¬ 
ified that this resolution is sufficient to obtain accurate 
results. We also verified that our simulation conserves 
energy with better than 1% accuracy. The thermal en¬ 
ergy lost by the crust is partially carried away by neutri¬ 
nos and partially conducted through the boundaries. 

Figure [5] shows the gradual evolution of the tempera¬ 
ture profile T(z) after the flare with the fiducial param¬ 
eters (see Section 2.2 and Figures 5, 6). During the first 
month, the initial peak of temperature at z ^ 500 m 
is reduced from ^ 5 x 10® K mainly due to neutrino 
losses. Then the peak continues to flatten and spread 
due to thermal conduction, forming a rather flat profile 
of T ^ 10® K in a few years. 

We find that plasmon decay and bremsstrahlung make 
the dominant contributions to neutrino cooling, and 
synchrotron neutrino emission becomes significant in 
stronger magnetic fields ~ 10^® G. Electron-positron 
annihilation dominates neutrino cooling only in the shal¬ 
low, low-density region of the crust, and its net contri¬ 
bution to the energy loss is negligible. 

With increasing Bz (and at fixed amplitude so of the 
waves excited in the flare) the deposited plastic heat 
increases, which increases the role of neutrino cooling. 
As a result, the relative contributions of plasmon decay, 
bremsstrahlung, and synchrotron neutrino emission de¬ 
pend on Bz- This dependence is shown in Figure[9l where 



Fig. 9.— Fraction of the post-fiare crustal heat lost through 
surface emission and various channels of neutrino emission, as a 
function of B = Bz . The flare is assumed to excite a pair of Alfven 
waves with so = 0.43 (Section 2.2) which are plastically damped 
in the crust. 
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Fig. 10.— Surface thermal flux caused by the plastic heating in 
the giant flare. Upper panel: fixed sq = 0.43 and varying B = Bz. 
Lower panel: fixed Bz = 3 x lO^'^ G and varying so- 



















































we also show the energy fraction that is conducted to the 
stellar surface and radiated away. The remaining energy 
fractioir (not shown in Figure 9) is conducted into the 
core of the neutron star. 

One can see from Figure |9] that only a small fraction 
of the stored crustal heat is conducted to and radiated 
from the stellar surface. For example, at = 3 x G 
less than 1% is conducted to the surface; roughly half of 
heat is lost neutrino emission and half is conducted to 
the core. 

The heat radiated from the surface produces a delayed 
afterglow emission of the flare. The solution of Equa¬ 
tion (1^^ gives the surface radiation flux F as a function 
of time. This flux is shown in Figure [TOl One can see 
that the surface flux peaks with a significant delay af¬ 
ter the flare — it takes the thermal conduction timescale 
(of months to years) to transport the crustal heat to the 
surface. 

The core remains much colder than the plastically 
heated crust. The heat conducted to the core can¬ 
not significantly boost its temperature because (1) the 
core has a large heat capacity, so a huge energy Fth 
1048^2 ijg requ ired to heat it to e.g. 10® K 

(| Yakovlev fc Pethickll200^ . and (2) the core is efficiently 
cooled by neutrino emission, and the cooling rate quickly 
grows at high temperatures. 


5. DISCUSSION 
5.1. Plastic damping and cooling 

In this paper we described the phenomenon of plastic 
damping of Alfven waves generated in magnetar flares. 
Our results may be summarized as follows. 

(1) Transmission. The flare generates magnetospheric 
Alfven waves with energy density 

Uo (23) 

where hb = is the tension of magnetic field lines, 

and So ~ 0.1 is the shear strain of the field lines. The 
waves are quickly transmitted into the crust of the neu¬ 
tron star. The transmission coefficient is T ~ 0.1 (Fig¬ 
ure 3), and most of the wave energy is transmitted after 
N ~ ~ 10 reflection events (Figure 6). The trans¬ 

mitted waves form a train of N oscillations propagating 
with velocity v ^ 10“^c and compressed by the factor of 
c/v (Figure 4). 

(2) Compression. The wave energy, which is initially 
spread in the magnetosphere, upon transmission becomes 
compressed. The energy density of the transmitted wave 
is 

U^-r-Uo, (24) 

V 

where v decreases to 10® cm s“^ as the wave train prop¬ 
agates toward the bottom of the crust. The transmission 
occurs at depths z of a few hundred meters where the 
crustal density p ^ 10^® — 10^® g cm“®. In this region, 
Uw exceeds the maximum energy that could be stored in 
the elastic deformation of the crust, Ue\ = ps\^l2, and 
the wave propagation is still sustained by the tension of 


magnetic field lines, ps- Therefore, the transmission also 
leads to the strain amplification: s^/s§ ~ U^/Uq ^ 10. 

(3) Plastic flow. The shear strain of the transmitted 
wave, s ^ {Tc/vfl^^so, exceeds the maximum possible 
strain of elastic deformation Sa ^ 0.1. Therefore, the 
wave induces a strong plastic flow of the crust, which 
dissipates the wave energy. A minimum plastic stress a 
is comparable to pscr and may be higher, as the shear¬ 
ing motion is very fast and the viscoplastic term rjs can 
dominate the stress. As the wave propagates into denser 
layers p ~ 10^^ g cm“®, the shear modulus increases to 
p Ri 10®®pi2 erg cm“® (Figure 1). The plastic heat den¬ 
sity deposited by the wave train is given by 

Uth asN ^ pScrsN. (25) 

The high Uth given by this estimate implies that the wave 
energy density Uw converts to Uth, i-e. efficient damping 
occurs. 

(4) Melting. Damping of the wave is buffered by melt¬ 
ing — plastic damping is inefficient where the heated 
crust becomes nearly liquid, and the wave continues to 
propagate to denser layers that have a higher rmeit(p)- 
As a result, a simple temperature profile T ra T^eiflp) is 
created by the plastic flow in an extended region of the 
crust (Figure 7). 

Most of the wave damping occurs at depth Zdamp 
where Tmeit is so high that the wave dissipation becomes 
marginally capable of melting the crust. Thus Zdamp is 
also the depth of the melted region. At this depth the 
following condition is satisfied, 

CvT^eit u^- (26) 

We found Zdamp ~ 500 m for a typical wave energy in 
magnetar giant flares and ~ 3 x 10^^ G (Figure 7). 
Deeper melting Zdamp 700 m is possible if the giant 
flare occurs in a flux rope of a particularly strong field 
B, > IQi® G. 

(5) Cooling. On a timescale of months to years, the 
deposited heat is mostly lost to neutrino emission and 
conducted into the core of the star (Figures 8 and 9). 
A modest energy Faft is conducted to the stellar surface 
and emitted in a delayed afterglow radiation. A typical 
energy radiated per unit area is Eait/S ^ 10®° erg cm“®. 
The timescale for the rise of afterglow luminosity is the 
thermal conduction time Fond ~ 10^ s. In a broad range 
of the flare parameters, the peak flux of surface afterglow 
is Fmax ~ (2 — 4) X 10®® erg cm“® s“^ (Figure 10). 

There are ways to refine our model of surface afterglow 
from plastic damping of magnetospheric waves. All sam¬ 
ple models shown in this paper assumed approximately 
vertical (radial) magnetic field in the upper crust. A 
strongly inclined field would significantly reduce thermal 
conductivity in the radial direction and delay the crustal 
cooling. It could also bolster a high crustal tempera¬ 
ture before the flare, which would give a deeper melted 
zone where plastic damping would be impossible. In this 
case, the flare could only cause heating of the deep crust 
where practically all heat is wasted to neutrino emission 
and inward conduction. Thus, a strong non-radial field 
component tends to reduce the expected afterglow emis¬ 
sion. 
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Our presented models assumed iron composition ev¬ 
erywhere in the crust, including the blanketing envelope. 
Light element compositio n of the envelope woul d increase 
its thermal conductivity (jPotekhin et al .11200,111 , decreas¬ 
ing the internal temperature and reducing the depth of 
the melted layer in the pre-flare crust. Therefore, if mag- 
netars have a light element envelope, their post-flare cool¬ 
ing occurs faster. This effect somewhat increases the af¬ 
terglow flux, especially at early times, and may offset the 
opposite effect of the non-radial magnetic field. 


5.2. Other mechanisms of Alfven wave damping 

Nonlinear interactions of Alfven waves in the mag¬ 
netosphere provide an a dditional damping mecha nism. 
The existing estimates (|Thompson k. Bla^ I1998D sug¬ 
gest that this mechanism will be dominant at very high 
amplitudes of the waves, so ~ 1- The nonlinear interac¬ 
tions occur as the Alfven waves bounce from the stellar 
surface and collide in the magnetosphere. The nonlinear 
terms in the electrodynamic equations show two types of 
wave interactions: 

(1) A-\- A ^ F-. two Alfven waves A convert into a fast 
magnetosonic wave F (which may escape the magneto¬ 
sphere) . The damping of Alfven waves by this “3-wave” 
interaction occurs on the timescale. 


■^/c ^ ( klY ^ 


(27) 


where A = 27r/fc|| is the wavelength along B (compara¬ 
ble to the length of the magnetospheric field line L), ^ 
is the characteristic displacement in the waves, and k± 
is the wavevector component perpendicular to the mag¬ 
netic field. The Alfven waves, which are ducted along 
the curved magnetic field lines, may be expected to have 
k± ^ fc||. 

(2) A + A A + A: two Alfven waves generate two new 
Alfven waves. This “4-wave” interaction initiates a cas¬ 
cade to high k I , which may lead to the w ave dissipation 
on small scales (jThompson fc Blaeslll998ti . The damping 
time due to this higher-order process is 


^ A/c 

{k^o^^ W) ■ 


(28) 


The time tdamp given by Equations (l?7l) and (051) should 
be compared with T~^Llc ~ lOL/c, the lifetime of the 
Alfven waves to transmission and plastic damping in the 
crust. The numerical coefficients in Equations (1271) and 
((551) have not been calculated, however the estimates sug¬ 
gest that if the flare generates sq ^ 1, the nonlinear wave 
interactions can reduce so to a value ^ 1 before the waves 
are damped plastically in the crust. 

The wave cannot be completely damped by the plastic 
mechanism. In particular, at strains |s| < Scr it prop¬ 
agates with no significant damping. The residual wave 
train will reach the bottom of the crust and enter the 
liquid core. It will travel through the core along the mag¬ 
netic field lines and after time ^ 2r/v (typically shorter 
than 1 s) the train will again emerge somewhere at the 
bottom of the crust and continue to propagate upward. 

The low-amplitude waves will continue to travel 
through the magnetosphere and the star for a while. 
Their lifetime at any given amplitude s is limited by the 


nonlinear interactions in the magnetosphere tdamp oc s“^. 
The Alfven waves are also subject to gradual ohmic 
dissipation, as their propagation involves excitation of 
electric currents demanded by V x B yt 0. After the 
flare, the effective resistivity of the magnetosphere is 
controlled by the threshold voltage of electron-positron 
discharge that organizes to condu ct the electric currents 
((Beloborodov fc ThomDSonll2007ll . 

5.3. Observed afterglow 

Sudden crustal heating followed by gradual crustal 
cooling was proposed to p ower the afterglow of t he gi¬ 
ant flare in SGR 1900-1-14 ((Lviibarskv et al.ll2002l l. The 
afterglow was extremely bright in the first hours after the 
flare, J§f ^ lO^’^—10^® ergs“^, and during t he next month 
it sh o wed a power law decay Af oc (iWoods et al.l 
l2001h . iLvubarskv et al.l ((200^ explored how heat should 
be deposited to give the observed afterglow light curve 
and found that heating should be approximately uniform 
throughout the 500-m-deep layer below the surface. This 
implies, in particular, enormous heating in the shallow 
layers z 100 m. The heating mechanism in the low 
density layers is unclear and certainly cannot be provided 
by plastic dissipation. Therefore we do not attempt to 
explain the early afterglow of SGR 1900-1-14 by crustal 
heating. We also note that the afterglow spectrum was 
nonthermal ((Woods et al.ll200ll) . which suggests a mag¬ 
netospheric source. 

Plastic damping of magnetospheric Alfven waves pro¬ 
duces a well defined temperature profile of the crust: 
T ft! T^e\t{z) down to Zdamp- This leads to specific pre¬ 
dictions for the afterglow light curves (Figure 10), with 
the surface flux F ~ (2 — 4) x 10^^ erg cm s“^ on 
a timescale ^ 100 d. This flux and timescale appear 
to be consistent with observations of some less energetic 
“transient” magnetars after their bursting activity. 

In particular, the luminosity of SGR 1627-41 after its 
outbursts in 1998 and 2008 showed a decay on a. year 
timescale (iMereghetti et al] 120061 lEsposito et al.l 120081 
lAn et all l2012f ). The luminosity at t ~ 100 d was 
.if ~ 7 X 10®'^(d/ll kpc)^ erg s“^ after the 1998 out¬ 
burst and .if ~ 2 X 10®'^(d/ll kpc)^ erg s“^ after the 
2008 outburst, where the distance d ft 11 kpc was in¬ 
ferred from the apparent location of SGR 1627-41 in a 
star-forming region (Hurley et al. 1999). The decay on a 
year timescale is consistent with the crust melting down 
to Zdamp ~ 300 m, and the observed luminosity .if is 
consistent with the melted crust area occupying ^ 10% 
of the stellar surface. 

Swift J1822.3-1606 provides another example. It pro- 
duced afterglow e mission following the out burst in 2011 
((Rea et al.l 120121 iScholz et al.l 12(112 . I2014( ). Similar to 
the afterglow of SGR 1627-41, its light curve may be 
described as a double exponential, with the second 
(longer) exponenti al component visible after ~ 100 d. 
IScholz et al.l ((2014H used a crustal cooling model to de¬ 
scribe both the early and late afterglow components in 
Swift J1822.3-1606. In their model, heat deposition is 
a phenomenological parameter adjusted to reproduced 
observations. We find that plastic damping of magne¬ 
tospheric waves is only capable of explaining the late 
afterglow component, and the early component must in¬ 
voke a different heat source. The late component has the 
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luminosity and decay timescale similar to those observed 
in SGR 1627-41, consistent with the crust melting down 
lO -^damp ^ 300 m. 

A reliable identification of the crustal afterglow is com¬ 
plicated by the presence of another, nonthermal, emis¬ 
sion component. The nonthermal sourc e is likely presen t 
during the afterglow of SGR 1627-41 (lAn et a.1.1 l2012h . 
and nonthermal hard X-rays are unambiguously de¬ 
tected in the transient magnetar IE 1547.0 -5408 during 
its a f terglow following th e 2009 outburst (jEnoto et al.l 
120121 : iKuioer et al.l 120121 1 . The nonthermal activity is 
usually associated with the twisted equilibrium mag- 
netosphere, which car r ies persistent elec tric currents 
(|Thompson et al.ll2002t l^eloborodo^ 1201311 . The twist 
is ohmically dissipated over a year timescale, which hap¬ 
pens to be comparable to the timescale of crustal cooling. 

Another complication is the expected external heat¬ 
ing of the stellar surface bombarded by magnetospheric 
particles. This heating occurs at the footprint of the 
current-carrying magnetic field lines (“j-bundle”). As 
the magnetosphere slowly unt wists, the i-bundle shrinks 
and so does its hot footprint (jBeloborodovI 1200^ . Such 
shrinking hot spots have been observed in several tran¬ 
sient magnetars, including the canonical transient mag- 
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netar XTE J1810-197. Eollowing an outburst in 2003 it 
showed an X-ray afterglow decaying on a year timescale, 
with luminosity .if ^ 2 x 10^^ erg s“^ at t ^ 1 yr 
(j Gotthelf fc Haloernl l2007| l. The observed area A(t) 
and luminosity of the hot spot evolved in agree¬ 

ment with the predictions of the untwisting magneto¬ 
sphere model. Similar shrinking hot spots were ob¬ 
served in IE 1547.0-5408, GXOU J164710.2-455216, 
SGR 0501 -f4516. SGR 0418+ 5729 (see the data col¬ 
lection in iBeloborodovI (j2011[l and refere nces therein) 
and more recently in Swift J1822.3-1606 (jScholz et al.l 
l2014fl and the Galactic C enter magnetar SGR J1745-2900 
l|Coti Zelati et al]l2015ll . 

Strong Alfven waves and deep plastic heating are cer¬ 
tainly expected in energetic events, in particular in giant 
flares. All three giant flares observed to date were emit¬ 
ted by persistently active magnetars, which maintain a 
high level of both magnetospheric activity and surface 
luminosity. It is possible that plastic damping of Alfven 
waves is the main mechanism that keeps the crust hot in 
these objects. 
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